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ABSTRACT 


Low sneed / low manoeuvrability aircraft are currently quite susceptible to being killed 
in attacks by the ubiquitous infrared mussiles. A theoretical analvsis applied to an en- 
counter simulation seems to indicate that it 1s possible to use the infrared jammer to 
defeat second generation infrared mussiles. The theoretical analysis of a simplified case 
of a conical scan reticle with frequency modulation jamming leads to expressions for the 
target’s position, as seen by the mussile seeker, under no-jammuing and under infinitely- 
powerful-jamming conditions. The intermediate-power case is dealt with by numerical 
analysis for a selected, non-optimal situation, as the closed form solution 1s not 1imme- 
diately apparent. The analvsis indicates successful jamming in the situation studied. In 
the scenario Where the infrared missile 1s an almost continuous threat during the air- 
craft’s flight, infrared jammers and low visual signature paints, and perhaps low infrared 
signature paints, are short-term solutions that are potentially useful in increasing the 


survivability of these aircraft by reducing their susceptibility to infrared missile kills. 


Cy it 


DISCLAIMER 


The reader 1s cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within the 
time available, to ensure that the programs are free of computational and logic errors, 
thev cannot be considered validated. Any application of these programs without addi- 


tional! verification 1S at the risk of the wsen 
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I. INTRODUCTION 


A. BACKGROUND 

When low speed / low manoeuvrability aircraft operate from bases which cannot 
be secured to a sufficiently great depth (more than several kilometres), they are vulner- 
able to surface-to-air missiles (SAM). As small, man-portable SAMs are quite readily 
available, they constitute a high probabihty threat. 

When out of the SAM engagement envelope, air-to-air missiles (AAM) remain as 


potential threats to these aircraft. 


B. SCOPE 

This thesis will address a particular aspect of this threat - the threat by infrared (IR) 
guided missiles. This is because infrared missiles are ubiquitous and rather successful, 
as can be seen in Table | on page 2 ((General Dynamics,1988), (Loral,1989) and 
(Nicholas,1988)) and Table 2 on page 3 (Loral,1988,Slide 87329).1 

For this thesis, we will consider aircraft that are already deployed (as opposed to 
an aircraft still to be designed). Therefore the scope of survivability enhancement by 
vulnerability reduction 1s more limited than increasing survivability by susceptibility re- 
duction.2 
Hence we will concentrate on the (hopefully) easier task of finding measures to reduce 
the susceptibility of these aircraft to being hit by an infrared mussile. 

Due to U.S. Government restrictions on dissenunation of security-sensitive infor- 
mation, this thesis is required to be unclassified. This constrains the thesis to a more 


generic treatment of the problem. 


C. GENERAL SCENARIO 

For the purposes of this thesis, the scenario is that the aircraft launches from a base 
that can onlv be secured to a depth of a few kilometers. 1.e.,the aircraft 1s within the 
acquisition and launch envelope of threat missiles (especially easily-concealed man- 


portavle SAMs) on takeoff as well as during the landing phase of flight. (An example 


1 Quoted from Moore D.,(OSD), 23rd [RIA, TRIS 1985 IRCM Symposium, p.1 A-2, 1985. 


2 The terms used here are from Ball (1985). The probability of kill of an aircraft ( P, ) 1s the 
product of probability of hit by a threat weapon ( P,) and the probability of kill given that the air- 
craft is hit ( P,,,). Susceptibility 1s associated with the probability of hit whereas vulnerability is 
probability of kill given a hit by the threat weapon. 


Table 1. QUANTITIES OF MISSILES 
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Table 2. KILLS OF AIRCRAFT BY WEAPON TYPE 
(Note the large fraction of infrared missiles.) 











of this aspect of the scenario is an Antonov An-26 or Antonov An-32 aircraft operating 
out of Kabul.) The aircraft then climbs to its cruising / operating altitude, typically be- 
tween 5 and 30 thousand feet. While it 1s at medium altitude, the more likely infrared 
mussile threat is that of the AAM. 

Additionally, it is assumed that the launch platform of the infrared missile can reach 
a position such that the aircraft is within the acquisition and or launch envelope of the 
missile. The aircraft, which is usually unarmed, is unassisted by other aircraft and lacks 
the speed to run away from an airborne threat. 

Since infrared missiles are fire-and-forget weapons, multiple infrared missiles can 
attack the same target. Therefore, this threat must also be handled. Hlowever, we will 


also assume that the multiple infrared missiles are of one tvpe, during one attack. 


I]. THE AIRCRAFT CHARACTERISTICS 


A characterization of the aircraft in terms of the factors significant to the infrared 


missile countermeasure problem are as follows:- 


A. GENERAL INFRARED SIGNATURE FACTORS 

The infrared signature of an aircraft 1s dependent, to various degrees, on a number 
of factors. Various energy sources (and sinks) cause heating (and cooling) of the aircraft, 
and also reflect off aircraft surfaces. The energy is distributed over the aircraft primarily 
by conduction and radiation (and convection, in the case of the engine exhaust plume). 
The surface temperature and the surface spectral emissivity together determine the 
emitted power. The incident energy, geometry and spectral reflectance determine the re- 
flected energy component of the infrared signature. A summary list of the various coni- 
ponents follows (Higby.1972). (Where available, typical values are also quoted 
(Loral.1988,Slide 91537).) 

1. Emission component energy sources/sinks 


e Engine (thermodynamic cycle) 
(11ot metal typically evits 1] KW/sr/um at > 1.5 xm) 
(Engine exhaust shroud i. pically emits 100 W/st/um > 3 um) 


e Engine exhaust plume (Tvpically emits 50 W/st/um at 3.8 xm to 4.8 wm) 
e Avionics/Heat pumps. 


e Solar heating 


e Acrodynamic heating via adiabatic effects and friction ’viscous effects. Aerodvna- 
mic heating Is not significant at the speed at which these aircraft fly. ( < 4O0Okts 
TAS). 

e Radiation exchange with atmosphere and via nocturnal cooling. 

¢ Jree stream conduction (Skin surface emission 1s typically 10 W/st/uin at 8 um to 


14 ym) 


2. Surface emissivity 
Total emissivity (as opposed to spectral emissivity) is around 0.9 for the more normal 
painted surfaces. 
3. Reflection component energy sources 
In addition to the emission component energy sources (which are also reflection 


component energy sources), the other energy sources include:- 


e Sun - Sunhght reflecting off aircraft surfaces, particularly when specular reflection 
occurs, 1s known as ghnt. At typical values of 100 W fst/um at <3 um, sunhght 
reflection 1s a significant contributor to the infrared signature, when it 1s present. 


e Skv - This is indirect sunlight, scattered by the atmospheric components (e.g. 
clouds). 


e Ground (earthshine) - Again this is indirect sunhght, 1.e,,sunhght reflected off the 
earth’s surface or infrared radiation coming from the heated earth’s surface. 
4. Reflectance 
These values are highly variable in the infrared bands for painted surfaces, such 
as these aircraft would have. 
5. Geometry/Structure 
Whuile not directly a source of energy, the geometry of the aircraft does deter- 
mine how the infrared energy 1s spatially distributed. Geometrical effects include:- 


¢ Multiple reflections - The “energy packets” reflect off several surfaces before they 
reach the infrared detector. This redirects the infrared energy in complex ways. 


Absorption/re-enussion sites - Where a site 1s in relation to the energy source will 
determine whether it will absorb energy by radiation and/or convection. 


e Shielding - The energy source is masked by aircraft structures between the source 
and the infrared detector. 


e Shadowing - Part of a surface that appears darker because the energy source is 
shielded from that surface by intervening aircraft structures. 


e Engine exhaust plume shape and size - This 1s a function of the engine operating 
conditions and the airflow field around the aircraft. 


¢ Aerodynamic heating - This arises from adiabatic compresssion of air and from 
drag. Where this occurs, and the degree to which it occurs, is again determined bv 
the shape, size and position of the aircraft structures. 


6. Atmospheric transmission 
The atmosphere is a selective absorber. (Cooper,!1988.p.1-6) In the infrared 
bands, there are two windows, one at 3 um to 5 wm and the other at 8 uni to 14 um 
waveband. Most of the detectors used in infrared nussiles seekers are designed to work 


in these windows. Most detectors are still using the 3 wm to 5 ym waveband window. 
“% 


B. INFRARED SIGNATURE VS. PHASES OF FLIGHT/MISSION 
The takeoff and chimbout infrared signatures of an aircraft are significantly greater 
than during the cruise, on-station, descent or landing phases. The exhaust gas temper- 


ature varies by a factor of 2 or more between these phases of flight. 


wi 


C. INFRARED SIGNATURE MODELS AND MEASUREMENTS 

Many infrared signature models exist. However, the more readily available ones 
only model the engine and only the jet-engmes eg PIREPS (GEC, 1976) and Asien 
(ASD-WPAFB, 1975). There was no trace in the open literature of infrared models 
dealing with the turboprop engine, which typically power these aircraft. A few, (e.g, 
HIDE (Fhgby, 1972)) do consider the other factors such as skin surface emissions. Nel- 
ther infrared signature models nor measured signatures of aircraft were available for 
inclusion in this thesis. It should be noted that there is a lot of similarity of these models 
to the surface rendering models found in CAD,;CAM applications. 

The true infrared signature of the particular aircraft, which we desire to protect, 
needs to be verified and or determined by actual measurement with appropriate 
radiometers. (In the analysis that follows, only the relative radiant intensity of the 
jammer with respect to the aircraft 1s required. However, the absolute power of the 


jammer can only be specified if the aircraft infrared signature values are known.) 


D. VULNERABILITY 

Since vulnerability of an aircraft is defined as the probability of kill, given a hit, the 
definition of “kill” strongly influences the value of P, . Aircraft kill would be when the 
aircraft is destroved. However, if a kill is defined as denial of mission, then the P, may 
be larger. since electronic equipment, other support equipment and human operators are 


all susceptible to serious mission-hindering damage. 


E. MANOEUVRABILITY 

Low speed / low manoeuvrability aircraft are quite far, in flight performance, from 
an agile fighter aircraft. The climb rate of the low speed / low manoeuvrability aircraft 
is typically of the order of a few thousand feet a minute and their speed 1s in the low 
hundreds of knots. Acceleration normal to the aircraft’s forward flight is limited to a few 


times gravity. This prohibits evasion of contact by acceleration and jinking. 


Il. THREATS AND THE ENVIRONMENT 


A. INFRARED MISSILE SYSTEM 
Here are some typical features of infrared missiles. 
1. Airframe and propulsion 
Usually the infrared missile has a cylindrical body with cruciform delta and/or 
trapezoidal wings and canards. The nose is blunt and transparent at the tfrared 
waveband of interest. Propulsion for AAMSs tends to be from single stage motors. SAMs, 
on the other hand, are often two-staged motors, with a booster to get the missile up to 
flving speed rapidly and then a sustainer motor for continuing the flight. The flight per- 
formance is characterized by the acceleration and velocity vectors attarnable as a func- 
tion of time, and the maximum ranges this gives. 
2. seeker 
a. Sensor 
Earlier generation detectors were made with lead sulphide. Indium 
antimonide is more widely used nowadays. Lead sulphide detectors work around the 3 
um waveband whereas the detectivity of cooled indium antimonide peaks at around 5 
um. This allows the cooled indium antimonide detectors to work with radiation from all 
aspects with respect to the aircraft. This 1s in contrast to lead sulphide which can only 
work with the hot engine parts and engine plume. (Cooper,1988,pp.2-22,7-21) 

To accomodate the wide range (> 70 dB) in the scene and target irradiance 
as the missile flies towards its target, the detector is normally backed by an amplifier 
with Automatic Gam Control (AGC), so that the downstream processors have a nor- 
malized signal to work with. Significant parameters here are the attack and decay re- 
sponse times of the AGC and the saturation level. 

b. Optics and stabilization 

There is usually a telescoye in front of the detector to gather the infrared 
energy. To maintain the seeker line-of-sight with narrow beamwidth optics, the telescope 
is normally decoupled from the missile body motion. 1.e., the telescope is stabilized in 
inertial space. Significant parameters here are the maximum slew rates and the gimbal 
limits and the instantaneous field of view which is a function of the focal length of the 


optics and the size of the detector. 


c. Target discrimination|position determination 

The determination of the target existence and position in the presence of 
background scene ‘noise’ 1s accomplished by various means. Early generation missiles 
had rotating reticles. In this scheme, information 1s extracted from the envelope of the 
detected signals. 1.e.,amphtude modulation (AM). This was followed by stationary 
reticles with the scene moved (nutated) over the reticle. Information, in this case, 1s 
carried in the frequency modulation (FM). Both these schemes transform spatial infor- 
mation into time and lose spatial discrimination. More recently, pseudo-imaging sensors 
have been used that scan the scene with a spot (or spots). For mechanical convenience, 
rosette scan patterns are used. The size of the spot is usually of the order of magnitude 
of the target image size. In the future, even more information extraction will be possible 
with imaging systems that can determine target features within the whole target. The 
effort here is to increase the ability of the sensor to determine smaller signature targets 
at longer ranges in the presence of higher scene ‘noise’, ( from both natural and inten- 
tional ‘noise’.) 

3. Guidance 

The guidance law for many (perhaps most) infrared missiles nowadays Is pro- 
portional navigation of some sort. The proportional navigation ‘constant’ may vary 
during the flight of the missile. 

The tracking loops are used to control the missile. The details of these deter- 
mine the missile response to aircraft manoeuvring and target position noise (e.g., glint). 
The usual tradeoff between controllability and stability apphies. 

4. Steering 

Steering 1s accomplished by one or a combination of aerodynamic surfaces or 
vectored thrust. To save mechanical complexity, some missiles can only apply normal 
forces in one axis relative to their body coordinates. and depend on the rolling of the 
body to align the force mm the direction desired. 

5. Fuze 

Most AAMs are both proximity fuzed and contact fuzed. The field of view of 
the fuze and the triggering range and velocity determine the detonation point. Radar and 
infrared (coherent and non-coherent) fuzes are emploved. 

Smaller man-portable infrared SAMs are contact fuzed. Man-portable missiles 
need to be light in weight. Putting in a more complex, and hence larger, fuze, would re- 


quire a reduction in Warhead size, for a given weight linut and missile maximum range 


(amount of propellant). The smaller warhead would then be less effective in killing an 
aircraft at the longer range at which the fuze could trigger. 
6. Warhead 
Missile warheads are usually high explosive blast or high explosive 
blast fragmentation warheads. The pattern and the nature of the damage mechanism 


G@ictate how the aircraft will fare in a near miss. 


B. MISSILE ENGAGEMENT PHASES 
The engagement process goes through the phases of detection, acquisition, launch, 


flyout and detonation. Intervention in the engagement process is possible at anv stage. 


C. LAUNCH PLATFORNIS 
The launch platforms for SAMs range froin a man to vehicles. They can be very 
difficult to locate. Usually AAMs are carried by fighter aircraft. though, of late, they 


are increasingly being carried by helicopters. 


D. OPERATIONS 
The method of operation of the missile, e.g.,where the threat 1s likely to be, when, 
etc., are significant, but the combinations are manv and thev depend heavily on 


intelligence-derived information. This topic will not be discussed further here. 


E. ENVIRONMENT 
The urban environment that is given in the scenario means that the background 
clutter 1s higher than it would be in a rural environment. This makes missile engagements 


of targets more difficult in such environments. 


IV. RESPONSE 


In the scenario of this thesis, the aircraft is highly susceptible to successful attack 
by infrared missiles. Countermeasures have therefore to be devised that can increase the 
likelihood of the aircraft continuing its mission. Some possible countermeasures are now 


discussed. 


A. LAUNCH SYSTEM DESTRUCTION 
Destruction of the infrared missile launcher 1s difficult because man-portable infra- 


red SA Mis are easily concealca 


B. LAUNCHER ACQUISITION ENVELOPE AVOIDANCE 
In the scenario given, avoidance of the acquisition envelope of the launcher is not 


possible due to the limited geographic depth assumed. 


C. LAUNCHER ACQUISITION ENVELOPE REDUCTION 
Reducing the acquisition envelope of the launcher is possible because the acquisi- 
tion 1s often visual for infrared missile engagements. Low visual signature paint schemes 


are required. 


D. THREAT PROPAGATOR DESTRUCTION 

Destruction of the missile in flight can be achieved either directly by directed energy 
weapons (DEW) or by triggering the fuze of proximity fuzed missiles. The DEW ap- 
proach involves first detecting the mussile, tracking it, and then attacking the airframe 
or trackingfguidance system, or the fuzing system. DEWs (laser, electromagnetic or 
charged particle beam) are not vet available for operational use. Electromagnetic pulse 
triggering of the warhead 1s also a future option. Pointing of those weapons is also 
problematic because the missile has to be tracked first. If the spot scanning seeker proves 
difficult to avoid (see below) , then it may be necessary to make this a viable option. 
Proximity fuzes tend to be sensitive in the direction perpendicular to the longitudinal 
axis of the missile so that they are difficult to access for self-protection i.e.,from the 
forward direction. While it is possible to trigger them, it may not be straightforward to 


put enough signal into the fuze early enough to be useful. 
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E. THREAT PROPAGATOR AVOIDANCE 

The threat propagator can be avoided by out accelerating the missile and/or by 
manipulating the seeker’s tracking ability. It 1s very unlikely that low speed / low 
manoeuvrability aircraft can out-manoeuvre a missile because these aircraft cannot 
generate a large enough acceleration vector and/or displacement from the missile’s flight 
path. Also, this method would require a reliable means of determining the presence and 
location of the missile. Manipulating the seeker can be achieved by reducing the signal 
provided to the seeker, increasing the noise at the seeker or providing other (false) tar- 
Bers for the seeker. 

I. Reduce signal 

This can be accomplished by reducing the aircraft’s infrared signature. This can 
be done internuttantly, e.g.,bv engine power reductions at critical points in time. or 
continuously. For the latter, hot parts can be shielded (insulated) and airflow altered to 
mix with and or cool the hot parts. This will require structural and aerodynamic work 
(© execute. 

For reflected energy reduction, low infrared signature paints can be considered. 
A potential problem with these paints is that the thermal load on the aircraft will in- 
crease (Ball.1988,p.293). 

Another technique is *v ‘capturing’ the Automatic Gain Control (AGC) of the 
amplifier of the mussile’s detector. Essentially, this involves depositing enough energv 
@methe detector for a brief period (of the order of the attack time of the AGC) to cause 
the gain of the detector amplifiers to decrease to the extent that the real target signal 1s 
lost in the noise. The period of energy deposition has to be brief. or else the real target 
signature 1s actually enhanced, which is the opposite of the intention. To keep the de- 
tector blind for a prolonged period, this energy pulse is repeated before the AGC can 
decay completely. The overall intent 1s to keep the signal available to the missile proc- 
essing electronics as low as possible. The technique works particularly against missile 
seekers that derive information from the envelope of the signal, as the first generation 
AM seekers do. 

2. Increase noise 

Noise can be increased either by tactical use of natural sowrces, such as bv 
placing the target up-sun of the missile launch platform, or by dehberate use of high- 
powered infrared sources. For the latter case. if the noise source does not saturate the 
detector, it may very Well provide a much stronger signal for the missile to home on, 1.€, 


become a beacon. 


3. False signals 

Flares are a commonly used means of countering infrared missiles. They are, 
however, expendables. In the case of these aircraft which will remain (because they 
generally have low excess power) for several minutes in the acquisition envelope of the 
infrared missile, reliable warning of the presence of the missile must be available to 
trigger flare release. Current false alarm rates and probability of detection of present 
day missile warning equipment are not adequate for this task. Because of the lack of 
manoeuvrability of these aircraft, separation of the flare from the aircraft, if needed, will 
not be aided greatly by manoeuvre. Rather, the ejection speed of the flare and gravity 
will be the determinant of the rate of shift of the centroid of the target and flare signa- 
tures.3 Moreover, flares with the combination of the right spectral quality (to defeat 
multi-colour seekers), quick burn build up and long burn time with low weight, are not 
readily available. Flares may also be difficult to use on takeoff, if the takeoff path is over 
populated areas. 

Modulation of an infrared source to provide misleading signals to the seeker 
can be done. (Devices used for this purpose will be called infrared jammers.) Against 
AM seekers. in addition to the AGC capture, the phase of the envelope variations can 
be manipulated by the infrared jammer and so the position vector determined by the 
seeker will be false. We do not deal with these seekers here as it 1s assumed that they are 
relatively easy to deal with. The second generation musstie seekers tend to be conscan 
(FM) types. An analvsis done to determine the potential of the infrared jammer to de- 
ceive such a seeker has been carried out and shows that it potentially can do so, as will 
be discussed later in Section “VY. ‘Conscan (FM) seeker vs. infrared jamm€r ferrin. 
thesis. Because of the narrow instantaneous field of view of a spot scanning seeker, it is 


unlkely that significant false signals can be inserted into these seekers in this way.4 


3 It may, on occasion, be useful to have a low separation speed against a spot sceker missile 
that uses kinematic 1ollowing to try to defeat flares. In that case, the missile miss distance will be 
small. Therefore, this can only work if the missile does not have a proximity fuze. 


4 While the spot secker missile must be dealt with, readily available information on these 
missiles is lacking. Therefore, with the time constraints of this study, we can only speculate on 
possible countermeasures against these missiles. 


V. CONSCAN (FM) SEEKER VS. INFRARED JAMMER 


A. BACKGROUND TO THE ANALYSIS 

This 18 a reticle in which the image is moved (nutated) over a stationary reticle. The 
resulting frequency modulation of the signal out of the detector can be processed to ex- 
tract the 2-D position of the target relauve to the reticle centre. 

The following analysis is derived from papers by Gedance(1961) and Suzuki(1979). 

Gedance’s analysis is for a target 1n the form of an illuminated disk of uniform and 
constant irradiance, and a spoked (or Wagon wheel) reticle with a transmissivity function 
that is either one or zero as a function of the azimuth of the reticle. His analvsis stops 
at the determination of the instantaneous frequency generated by the target image on 
feereticie, Me docs ot derive the 21D position of the targct. 

Suzuki considers a point target for a reticle similar to that analysed by Gedance. 
Suzuki then shows the relations necessary to obtain the 2D position of the target. 

Neither Gedance nor Suzuki deals with the jammung case which the following 


analvsis does include. 


B. ANALYSIS 

In the following analysis, we want to determine the apparent target position, as seen 
by the sceker. while under jamming. To do so, we obtain, in order, the flux seen bv the 
detector, then the instantaneous frequency of the flux signal and, finally, the expressions 
for the target position, as determined by the seeker. The expressions for the target posi- 
tion will be for the two extreme cases of no-jammer and for an infinitelyv-powerful- 
jammer. Onlv one case is dealt with here for the intermediate-jammer-power range. This 
is because the general solution to this case 1s not rmmediately apparent and hence a 


numerical analvsis had to be used. 


© FLUX SIGNAL 

Assume a uniformly illuminated target disk of radius 6 moving in a circular path 
of radius r about the origin Opo; (See Figure 1 on page [4.) The position of this origin, 
With respect to the reticle centre O,,; , 1s in fact the “position” of the target with respect 
to the sceker boresight. The target image centre 1s at (D, a) with respect to the reticle 


Genie One; - 
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Target Image 


OnoT 


Figure 1. Reticle-Target Geometry 


In general, the flux, /(v) , transmitted through a reticle is given by 


F(t) = 12) | 1(s, 0)\dA (1) 
A 


where 
I(t) = irradiance of the image, 
A = area of the image, 
(s, 0) = polar coordinates relative to the reticle centre Opry, 


T(s, 0) = transmissivity of the reticle. 
The equation of the circle that describes the target image is given by 


a ao cos(a — @) (2) 


Solving for s , in the above equation, we get the two roots 


s, = Dcos(a — 6) — ,/ D* cos*(a — 9) — (D* — 5°) (3) 


s) = Dcos(a — 6) + ,/ D’ cos*(a — 8) — (D’ — &°) (4) 


: a 
The elemental area, dA , is approximated by a rectangle of width ( eS ae and 


side Jengths (s, — s,) . Therefore, the incremental area 1s 


dA =(5)- 5 | ( -- Jeo | (5) 


maom Eqns. (1), (3), (4) and (5), the flux F(z) is now 


a+ aresin( pias ) 
D 


F(t) = (2D T(s, 8) cos(x — 6),/ D? cos*(a — 0) —(D? — 5*) d0—(6) 


1— aresin{ ccs ) 
D 


{f the reticle transmittance is independent of the radius, s , 1.e.,the reticle 1s a radially 


Saatnetric reticle, then 


I(s, @) = 1(8) (7) 


2m 


rue acomepeadescribedias.a Fourier 


And if 7(@) is a periodic pattern of period 22 (= 


Series thus: 


OO 


Cy = = - ) 6 cos( nat + én) 


a (8) 


P Go 
= = abe » fr cos(im0 + d,) 


n=! 





where 


s 
m 


= T(A)e 2" a0 (9) 


Sa 


Inserting this description of the reticle transmittance function (Eqns. (7), (8) and (9)) in 
the flux expression (Eqn. (6)), and reversing the order of summation and integration, 


We can write 


FU) = Ko} bot Y Ag(o (10) 
n=] 
wahene 
a+ arcsin( < ) C — 
Ay = 2D ae cos(« — 8),/ D* cosa ~ @) ~ (D* — 5°) dé (11) 
= arcsin( < ) 
and 
; 1+ arcsin( + a ae 
A,,(0) = 2D | ¢, cos(ain? + ,) cos(a = @).7 D” cos (a — 6) —(D" = 0 \ia a 


~a— arcsin ( - 
{— — 
dD 


Evaluating A, (Eqn. (11)), we get 


A+ arcsin( & (a5 oe =. Oe 
Ane cos(a — @),/ D® cos (a — @) — (D° — 0°) dé (1:3) 
a— arcsin = 
Let 0'=a—-—86 . Then get 
- arcsin( — ) a5 =e Eee 5 z xs 
= 1a cos(6’),/ D* cos’ 6’ — (D° — 5°) d{ — 8’) (14) 
arcsin =) 


Since the integrand is even, 


aresin( O_ ) 
D 


Ag = 2D6q, cos(0") © y cos’ — Ga) + 1 dé’ 
(15) 
arcsin( > ) : 
= 2D0dcp cos(@’) j/1— €. y sin’0’ d@’ 
0 . 
From tables (Beyer,!987), Eqn. (15) becomes 
: ’ ney ae ee : : won 
Ag = Dic] SAE J - € ) sin’O’ + a os aresin( = sin o) fro D 
(16) 
= 28 
Wehine 
A “Oe? 
A= yj =e (17) 
Likewise, for A, (Eqn. (12)). 
1+ aresin( 2 ) 2 a eee 
A,(t) = 2Dc 4 cos(nnO + $,,) cos(a — 8)! D* cos*(x — 8) —(D* — 6°) dO (18) 


I= arcsin( a ) 
D 


For simplification. change the integration variable to @’=a-—8@ . Also, expand the first 


cosine term in Eqn. (18). Then 


- arcsin( <. ) [Pyrcote (PFs ya 
emer 21), cos(7(a — 8’) + d,) cos B oN eM cos’ 6’ ‘e +1)d 
“arcsin( =) 
| —— ee ae 

= 2D6c,, cos(nnx + @,) cos nm’ cos afi _ .- y sin? @’ 


= arcsin( OL ) 
D 


+ sin(vnz -+- d,) sin 278" cos afi — (2 i sin’ 0’ jo 


arcsin ( 


(19) 


O 


The first term in the integrand of Eqn. (19) 1s an even function, while the second term 
is an odd function in @’ and the integral of an odd function taken symmetrically about 


the origin vanishes. Therefore, 


: 6 = ee 
; Ares = 
A,(Q) = 2Do6c,, cos(nmea + é,)| ( p ), cos @’ cos nin’ | — (2 : sin?6’ do’ (20) 
0 


It can be shown (Gedance.1961,Eqn. 9 and 10)5 that if <1] and we ] , then 

nm arcsin( -F ) gO). 2 nnd 

a 2cos 6’ cosnm@’ {1 — ? sin?@’ de'=J, (a) 

‘ 0 oe 
0 
where ia 1s the first order Bessel function of the argument C. 
Then 
= as Doc, cos(mmna + d,)J ‘Coe (22) 


From the geometry of the situation (See Figure | on page 14.), the triangle formed by 
the centres of the reticle (O,;;}. target and rotation of the target image about the reticle 


(Ore) ties 


2 ee 
Di = Ry+ = 2K cos(w,f =a 07) (233 
and 
r sin w,t — Rysin é7 
OF ON CECT heme eer ets ame (24) 
rcosw,t— Rycos 6+ 


& 
Establishing p as a normalized offset distance, —— , then 


D=r,/1—2p cos(wt —07) + p* (25) 
and 
sin w,— p sin 67 
v= Gnc0n| —— (26) 
cos w,l — pcos @r 


5 See Appendix A for a demonstration, by series expansion, of this approximation. 


A, can then be written in the condensed form 





A,()=A4,, cos O,(2) (27) 
where 
en / 2 nine 
Ay, = wee SX ]—- 2p cos(wr —0,) a 78 6c), —————— (28) 
r./1— 2p cos(wt -87) + pr 
and 


(29) 


A sin w,f— p sin 67 
eG — van arctan ——_—__—_——— ] + Py 


COS w,f— pcos 07, 


Bore that for p<l , A(7)= constant. Note also from Eqns. (27) and (29) that the infor- 

mation on the position (p , @;) 1s carried in the angle. i.e. the target image motion and 

reticle interaction gives rise to a frequency modulation (FM). If the reticle pattern is 
) Ue 


ick :) 
such that the conditions a <1 and ore = cain Deameto = then 


N 
FU)= MO Ag + > Anl2) cos O,(0 (30) 


n=1 
Where .\ is large but not violating the condition. 


D. INSTANTANEOUS FREQUENCY 

We now develop an expression for the instantaneous frequency of the flux signal 
as a step towards determining the information carried in the FM. 

The flux signal, F(z) 1s then bandpass filtered. Let’s call the bandpass filtered signal 
F(a). In general, a bandpass filtered signal can be represented in quadrature form or 


Semlpiex envelope form. 1.¢ 


? 


P(t) = P(t) cos wt — O(2) sin wt (31) 


6 A seeker with a wagon wheel reticle would normally be designed to have the target image 
of approximately the same size as the width of the transparent spokes of the reticle, at the radius 
at Which the target mage rotates on the reticle. This would give the greatest target-to-background 
discnmination. Also, dunng a subsequent bandpass filtering operation (which will be describe later), 
the higher order components (associated with 7 1n one of the conditions) of the reticle-generated 
sicnal are removed. Thus, in practice, in the vicinity of the design point of the reticle seeker, the two 


a 0 M210 
conditions of — <] and ——— 


< | Gan Usually be mmicl. 
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Or 
F(t) = Relg(t)e?"4 Ce 


Where, 1n this instance, if 771s the number of cycles of transmittance changes that a point 
on the image goes through in one rotation about the reticle, then w,= mw, . The two 


forms given in Eqns. (31) and (32) are related by 
Ee) tia eo) neD 


From the quadrature components P(r) and Q(z), we get the instantaneous frequency 
of that signal? (Roberts,1977) as 


PQQ — Q(2) PCa) 


- 34 
P+ O°(0 - 


cw {1) = as 


In a step reminiscent of synchronous or coherent demodulation of an FM signal, but 
using the scene rotation frequency as the reference frequency, we can extract the in- 
phase and quadrature-phase components of the signal. Later, we will show that this in 
fact gives the position of the target. To get an easily recognizable expression for the 
position in the end, we also include a scaling factor (=) , at this stage. The in-phase 


component, (t) , and quadrature phase component, g(t) , are described bv 


OXt)4+2 


i(t) ee: Frc. 210) een fa (35) 


Gt) J A Jers)ma sin wf 
w, 


Note that we can only find the position vector a half-cycle behind current time because 


the integration takes place over a full cycle. 


E. SIMPLIFIED RETICLE AND JAMMING SIGNAL 
It is rather tedious to continue the analvsis generally. So that we do not lose sight 
of the forest for the trees, we shall consider a special case of a reticle and of a jamming 


signal. 


7 Essentially the instantaneous frequency is given by the differentiation of the phase of F,(0) , 
with respect to tume, and offset by the carrier frequency. 


We can simphfy the reticle transmission function to a single sinusoidal azimuthal 
Variation, reducing the Fourier summation description of the reticle to a single term, 


besides the constant term. 1.€., 
(0) = + > cosl (6 — ,)] (36) 


In terms of the quantities defined in Eqn. (8), 


Cy = | 
l 
mi > 
Pp =— mg, 


moO iO tt = | 


iinen, irom Eqns. (17), (28) and (29), we get 
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Ape =, (37) 

A Oa | 2 10 i 
(0) a rr e Ory! l= 2p cos(w! —6-) sine J, ee (38) 

Aa — 2p cos(wi —87) + p 
and 

sin w,f — p sin @y ; 

SLi GOL) a |), (39) 
cos w,t — p cos 67 


We now have a description of the component of the signal that results from the inter- 


Pemon of the target and the reticle. in the seeker. 


F. VARIABLE IRRADIANCE FROM JAMMING 

To investigate the potential for jamming, we apply an optical signal combining a con- 
stant target irradiance with a sinusoidul jamming signal. The irradiance , /(t) , 1s then 
described by 


Kt) = Ie, + GLI + cos O2)] (40) 


where /,, 1s the target irradiance 


tgt 


and JL 1 + cos ©,(t) ] 1s the jammer irradiance which oscillates between 0 and 2/,. 


AN| 


Rearranging Eqn. (40) , we get 
{= 174 eos © 7) (41) 
Avinere 
1 = lige + (42) 


is the unvarying component of the irradiance. 


The flux signal for this special case 1s 
Fi) =, + J, cos O17) IL Ap + A, (2) cos ©1,(1) J (43) 


After high-pass filtering, the D.C. term, /,A, , will be rejected by the filter, and we get the 


bandpass signa] 


4 
Ft) = > Fy cos $4(0) (44) 
i 


Werke 
F, = [A,(1) P; = ©),(1) 
r= L Ag ¢.= ©,(7) 
Fy=1A()  6;= Of1) - O,(d) 
1s = PA) P, = © (1) a O,(1) 
In complex envelope form, 
F(t) = Relg(1)e/e| (45) 
ove de 
4 
gli) =) Fyel@no- (46) 


pe 


For compactness in subsequent expressions, let 


P(N 2 O(N) — Ot. (47) 


jap 


To find the instantaneous frequency of the bandpass signal, we need first to find the 
quadrature components (P(z) and Q(r)) of the bandpass signal. This is done using Eqns. 
(31) and (32) with Eqns. (46) and (47). We then get 


P(t) = Re{g(y)} 


4 
= Ne. cos D,(2) (48) 
=] 
and 
O(1) = Imig(2)} 
4 
= > Fr sin @,(2) (49) 
n=1 
Also, 
4 4 
P+ QW =) +d FyF cos,(1) — &{0) (50) 
k= (=1 


Taking the derivatives of the quadrature components described by Eqns. (48) and (49) , 


we get 
4 
Pi) =- » OF, sin D,(2) (51) 
n=) 
and 
4 
O(1) = » oor, cos D,(2) (37) 
n=] 
“b 
Also, combining Egns. (48). (49), (51), (52), gives 
4 4 
PUNO) — OWP() = >) ODF, F; coll) — &{0) (53) 


ei 


Substituting Eqns. (53) and (50) in Eqn. (34) , the instantaneous frequency in this case 


is described by 


4 4 
>» >, ODF F cos(,(0) — 00) 
«) (2) => G) > + a ee (54) 


4 4 
YY Fak; cos(@,(2) — (0) 


Now we define a jam-to-signal-plus-jam ratio, £,, and ee restrict the analvsis to the 


If p<1 (i.e., small 


} J 


7a as 


tet 





case of small target offset from centre. Let f= 


target offset from centre), then 
Ayp= A (OA, (SS) 


Hence 








E ma orem ae (S6) 


and we then get 


4 4 
> >, bf, col) — 040) 


joe Ihe 


w(t) = 0, += (Sa 
> >, Behr cos(Py(2) = (0) 
eit 4 
alnlend 
fl ff k= 
Pe) if kel (58) 
and 
| if l=1 
p -| | oF 
De eal 


G. SPECIAL CASES 

Let us now examine tlus result in the two limits of f£, =0, no jamming, and f,= 1, 
the infinite jamming case. 

1. No-jamming case 


Consider the no-jamming case when /,=0 . 1e.,when f,=0. Then 


i P 
= 0 if k#] (60) 
and 
ae 
iat if [#1 cy 


Sostitutine Eqns. (60) and (61) in Eqn. (57), wesget 


d,(1\(1 iP 


(1) 
- a, + b(t) — @, 


w(t) = ow. + 


oe (62) 
= Ul!) 


| —pcos(w,t — @7) 
SS) a a a ea 
1 — 2p cos(w,t — 87) + p 
Now that we have the instantaneous frequency of the signal, we will show that the co- 


herent demodulation of Eqn. (35) does mndeed give the position of the target. For the 


in-phase component, 


Or4+z 
O, qd 
= aly WANr COS w,tdl 
ae 
Ww 


id 


(63) 





Op +7 
] @, l—pcos(w,t — 07) 
me, (0) ae ee ean 1) 0) > w tal 


Or—7 ote (ne) ae 


Ww, 


Defining the change of variable @’ = w,: — 6, and expanding the trigonometric terms, we 


convert the integral into the form, 


ty 
A 


- 


cos 0’ cos @-— sin @’ sin @7;— pcos 67-cos "6" +p sin 8,sin @’ cos 6’ 


~~ 
| 
=i] 


- de’ 
1] — 2p cos(w,1 — 67) + p° 
re P (64) 
2r cos Or cos 0’ — (p'2) cos 20’ + (p‘2) 16" 
7U 


— 
— 


1 — 2p cos(w,t — 07) + o 
From the tables of integrals (Gradshteyn 1980) we can evaluate this definite integral, as 


_ 2rcos és 7 ; 
a a eS 

2 
rp COS O¢ (65) 


Rycos @r 


Which is in fact the x component of the target position. Similarly, from the 


quadrature-phase component, Moe & get 


which 1s the y component of the target position. 

Hence the target location in two dimensions 1s now determined. (Note that 
whether or not there is jamming. the missile seeker, not knowing anv better, will still 
perform the same procedure, described above, to try to extract the target position.) 

2. Infinitely powerful jammer 

Let us now consider the case in Which we have an infinitely powerful jammer. 
ic. f, = 1 This case gives an instantaneous frequency expression that 1s not simple. To 
gain further insight, we will restrict the bandpass filter to an even narrower band about 
the carrier frequency ( w, ) of the flux-derived signal. The upper limit of the bandpass 1s 
just above twice the carrier frequency. Mathematically, we make @, and @, zero. Then 


Ean. (37) becomes 


D,(2) + (1) + [,() + (2) |] cosl®, (1) — ©,(2)] 


W(t) = w, 2+ 2cosL®,(1) — 0,(1) J 


(67) 
= > [O,() ai O() | 


And if the jamming waveform is similar in form to the reticle modulation waveform (as 
given in Eqn. (39)), and the jammer’s equivalent of the reticle rotation frequency is 
w,, = w,, then the position coordinates of the apparent target location are : 

> (Rrcos 6-+ R, cos 6,) 


(68) 


N|— N|- 


(Ry sin 67+ R; sin 6;) 


where A, & @, are the jamsming parameters corresponding to R; & @; of the reticle. Note 
igemenicctive pointing direction of the seeker to the apparent target is that for the 


centroid of the combination of the target signal and the jamming signal. 


H. SOME OBSERVATIONS ON THE ANALYSIS RESULTS 

The jamming signal cannot dominate the target signal because the jamming signal 
is itself modulated by the relative motion of the jamming spot and the reticle. 

The analysis here 1s of a reticle that has a transmissivity function that varies grad- 
ually in a sinusoidal manner in the azimuth direction. For manufacturing ease, most of 
these reticles in real missiles will have more sharplv changing transmissivities. Howe'ver, 
to linut noise, the signal from the detector will eventually have to be bandlimited. Thus 
this analysis is a good first order approximation for the real missile. 

It should be noted from the analvsis that the conscan reticle seeker does not derive 
its information from the amplitude envelope. This being the case, AGC jamming tech- 


niques may not prove to be particularly effective against the conscan reticle seeker. 


I. SIMULATION 
Numerical analysis was performed to examine the intermediate cases of f, . For 
simplicity and computation speed, only the narrow band version of the analvsis was in- 


corporated in this simulation program. 1.e., 


> >, bf 6; cos(@,(2) — &{0) 
k=] /=1 
2 


» 


c{1) = WO, 


= 


BB, cos(D,(1) — Pz) (69) 


> 
~, 
i 


] 


D,(0) + 6; P,(1) + [d,(2) + @,(1) |B; cosl (2) — (1) J 
P= p: 0, cos{_, (2) — @,(2) ] 


= 10), + 


Eqn. (69) was then inserted into Eqn. (63) and the corresponding equation for the v-axis, 
and was transformed into computer code. At the same time, to incorporate the missile 
aerodynamics, the 2D position vector computation from the seeker was inserted into a 
version of the TACTICS IV program (Meisberger,!983) which was ported to a micro- 
computer. TACTICS IV is a time step simulation of a missile flight. It computes the 
relative positions and velocities of a target and a missile, takes the seeker-derived posi- 
tion (with noise added, if desired), applies a proportional navigation guidance law, and 
integrates the effects of a commanded acceleration to obtain the new position and ve- 
locity of the missile. The modified TACTICS IV, which we call TACTICS IV (IR) in- 
corporates the infrared seeker behaviour with an infrared jammer. The seeker routine 
was patched so that it will branch to the new infrared seeker routine when desired.8 The 
paramcters required for the seeker are essentially the parameters of Eqn. (69). In addi- 
tion, the tabulated output file of the missile and target flight profiles was modihed 
slightly to be more easily read by a graphics program.9 The listing of the infrared secker 
portion of that program can be found in Appendix B. 

An example of the result of an air-to-air missile flight simulation 1s in Figure 2 on 
page 29 . In this example, the missile is launched at the target from a position directly 
behind the target. The target does not evade, but continues in a straight and level flight 
path. Figure 3 on page 30 is the input data to the TACTICS IV (IR) program, for this 
example. In this particular case, the missile departs from the target in a spiral-hke path. 
Note the increase in muss distance to 700 feet with the infrared jammer on, as compared 
to the 3 fect case when there 1s no jammer.!0 

This studv should be continued with sensitivity analysis studies and optimizations 


to identify the optimal jamming parameters. 


8 This is done by making the ISEEK parameter of the input data table equal to 5. 


9 The postprocesser described by Meisberger (1983) was not ported to the microcomputer. 
Instead, a commercial package. STATGRAPHICS, was used for the graphical output. 


10 The no-jammer case was sumulated in a separate run, with £, effectively zero. 
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: Missile and target flying from 


Missile flight profile (with jamming) 


Figure 2. 


left to right. 
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0.20 0.80 1.50 2.0 Caso 2-3/7 2.75 4.60 0 

0.235 0.240 1.05 0.910 0.830 0.745 0.630 0.580 0 

1.04 1.04 1.04 0-93 0.86 0.90 0.90 0.87 0 

0.116 0.127 0.198 0.162 0.141 0.113 0.070 0.051 0 

0.0 4.0 6.0 8.0 10.0 12.0 14.0 30.0 50.0 

4.0 4.0 4.0 4.0 4.0 4.0 4.0 4.0 4.0 

TEST STT 

. 9 001 0.0 0.0 25000.0 0.9 0.0 
0 006 28.6 10000.0 0.0 25000.0 0.20 
0 011 0.0 0.0 1.0 6.00 0.00 
0 016 19Tac -10.00000 -0.000000 0.0 0.0 
0 021 0.1 0.0 0.0 0.1 0.0 
0 026 0.01 0.01 0.0 0.70 30.0 
0 031 0.0 0.0 0.136000 0.0 =15 
0 036 us) 15000 .0 0.001 0.0 aS 0. 
0 041 0.0 0.50 0.05 30.0 0.0 
0 046 6.0 0 0.0 0.0 20000.0 
0 051 6788.0 0.0 5.0 0.0000000005 100.0 
0 056 0.90 15.0 2300.0 0.00 21.8 
0 061 4.3 0.0 3.0 5590) Fes 
0 066 1.2 8.0 2.4 0.02 160.0 
0 071 2000.0 0.0 1.0 0.05 0.1 
0 076 0.0 0.0 0.0 0.0 0.0 
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Figure 3. TACTICS TV (IR) Input data for missile flight (with jamming) 
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VI. RECOMMENDED COUNTERMEASURES 


In summary, given the scenario considered here, the most promising short term 
countermeasure solution 1s the infrared jammer. The infrared jammer has a good chance 


of working if it has the following characteristics:- 


e For all aspect coverage, the infrared energy sources should radiate all around the 
aircraft. This would implv more than one source element, as there 1s no singie point 
on the aircraft that has equal visibility in all directions. 


e The power output of the jammer source, in the band that the missile detector is 
using. needs to be at least equal to the target signature in that band and in the same 
direction. 


e The jammer should be able to modulate its infrared power emitted such that the 
waveform is similar to what the seeker forms when it views a steadily radiating 
target. 

The analysis done here is only a paper study, using what are intended to be (hope- 
fully) typical parameter values. The actual signatures of these aircraft under various 
conditions must be measured to provide accurate input. The infrared jammer source and 
modulation scheme must actually be verified by test and evaluation with real hardware 
io deternune if the assumptions and analysis are correct. e.g.,can we find IR sources that 
can provide enough in-band energy on the seeker’s detector. We have not determined 
vet the optimal modulation scheme against the conscan (FM) seeker. Neither have we 
considered the optimal modulation scheme when there are more than one type of missile. 
Evaluation must be extended to these cases to ensure that modulation schemes, while 
effective against a particular type of nussile seeker, do not decrease the miss-distance for 
@umer seeker types. 

This study must also be extended to the third generation scanning spot seeker 


threat, Which 1s increasing with time. 


ot 


VH. CONCLUSION 


Low speed / low manoeuvrability aircraft are currently quite susceptible to being killed 
in attacks by the ubiquitous infrared missiles. A theoretical analysis applied to an en- 
counter simulation seems to indicate that it 1s possible to use the infrared jammer to 
defeat second generation infrared missiles. The theoretical analysis of a simplified case 
of a conical scan reticle with frequency modulation jamming leads to expressions for the 
target’s position, as seen by the missile seeker, under no-jamming and under infinitely- 
powerful-jammiung conditions. The intermediate-power case is dealt with by numerical 
analysis for a selected, non-optimal situation, as the closed form solution 1s not 1mme- 
diately apparent. The analysis indicates successful jamming in the situation studied. In 
the scenario where the infrared missile is an almost continuous threat during the aur- 
craft’s flight, infrared jammers and low visual signature paints, and perhaps low infrared 
signature paints, are short-term solutions that are potentially useful in increasing the 


survivability of these aircraft by reducing their susceptibility to infrared missile kills. 
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APPENDIX A. VALIDATION OF THE APPROXIMATION 


a 5 
We want to demonstrate that if <1 and 7 


arcsin 
an | ( ), cos 6’ cosnmé@’ f 1 — es ) sin?@’ a'=( =) (70) 
0 


Por convenience, let 


= |. then 


do 


al ose (71) 


and call the left hand side of Eqn. (70) , f(x). Now the first order Bessel function in 


Peles [Orin 1S 
~ 1% 24 
I(x) ==> een : (72) 
2 PRK +B)! 


eee eee Sg (73) 


Expanding the integrand 1n the left hand side of Eqn. (70), using the Taylor series, and 


integrating, we get 


EOS ioe a 1S] g — 
SS Sr S407 280n ~ 


elo.) ee) Os | RO oes 
Bee te eS)’. 


ls a )« , the terms in the second row of Eqn. (74) are small in relation to those in 


the first row and hence can be neglected. Computing the coefficients of x , 











(74) 








1 ee 3 eons 
Tee) = 797 *~ 746% ta93% 7* (75) 


~ 
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‘Comparing Eqn. (73) and Eqn. (75), if rune. =aal 


Sinem 
Soos(*) =I (x) | (76) 


which is what we set out to show. 
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APPENDIX B. TACTICS IV (IR) PARTIAL LISTING 


(CO 
Proeram listings of the infrared seeker 
Portion of TACTICS 1V CIR). 

Requires TurboPascal 5.0 and the 
TurboPascal Numerical Toolbox. 
MH HNN ANNAN RN NR) 
unit tac8s; 


interface 
uses tac8s; 


{$I floatdef.inc} 
procedure csj_seeker(time : float; 
var seeker_theta_prev, 
seeker_phi_prev,. 
seeker_theta_cur,. 
seeker_phi_cur : float):° 


procedure GetOutputFile_Namel(var OutFile : texts 
var FileName : string: 
default_filename : string); 


procedure sget_root_name_of_files(var FileName : strings: 
default _filename : string): 


implementation 


{$I-} { Disable I/O error trapping } 
{$R+} { Enable range checking } 


uses 
Intesrat. Dose Crt+ Common. Funct: 


thetac_j_default 
del_theta_default 
beta_j_default = 
foplteecions) = 1).0's 
max_move_step = 0.1; 
filename_csj_default = "CSJ.OUT'; 
focal_lensgtn = 57.0; 


Gq 
thetac T_init = 
1 
U 


4159265/33 
Ors 


var 
wire wir}. del_w_c» del_wir. beta_j : float 
x_tjJ» vyitj- range_tot. theta_tgt : float: 
x_tl» vitl. ranse_tstl. theta_tstl : float; 
war t cyc. t_cyc_j : float: 
mer ts del t 2 : float: 
u_lim. 1_lim float; 
del_theta : float; 
Gigms. £NO J start ¢ float: 
del_wir_rev : float: 
1 : integer; 
consl : text; 


rho- rho_step : float; 

tietrac 1. thetac_j : float: 
temp_beta : float; 

new_x, new_y : float: 

del_ranse;, del_«. del_y : float; 
FileName ;: string; 


var 
LowerLimit. UpperLimit : float; { Limits of integration } 
Tolerance : float; { Tolerance in the answer } 
MaxIntervals : integer: { Maximum number of subintervals used } 
{ to approximate the integral } 
Integral : float: { Value of the intesral ) 
Numintervals :; inteser; { Number of subintervals used } 
{ to approximate intesral } 
Error : byte; { Flass if something went wrong } 
{SF+) 


function betal(t 1: float) : float; 

var tt : float: : 

begin & 
beta := beta_j; 

end: {function beta} 


function theta_ol(t : float) ;: float; 

var numere denom : float; 

begin 
numer :* sin(worst)-(rho*sin(thetac_T) 
denom := cos(w_rt)-(rho*cos(thetac_T) 
theta_o :2 m*(arctan(numer/denom)-phic_r); 

end; 


is 
): 


function theta_j(t : float) . float; 
begin 


theta_j := m_isCarctan((sin€(w_oredel_weor)t)-(rho_jxsin(thetac _4)))/ 
(cos((w_redel_w_r)st)-—(rho_j™cos(thetac_4))))- 
end; 


function thetedot_o(t : float) : float; 
begin 
thetadot_o := m¥wor*"(Cl-rho™cos(w_rst-thetac_T))/ 
(l-2*rho*cos(w_rst=thetac_T+erhoxrho)); 
end; 


function thetadot_j(t : float) : float; 
var numere denom : float; 


begin 
numer := (l-rho_d¥cos((w_oredel_wor)™t-thetac_j)); 
denom := (l-2¥rho_j*cos((w_orxt+del_w_r d-thetac_j)¢*¢rho_isrho_j)3 
thetadot_j := m_isn(woredel_ wor )*(numer/denom); 

end; 


function ohi_gCt + Flost) ¢ float; 
var tj : float: 
begin 
t41 s= theta_i(t); 
phi_g : 8 ti-(myworst); 
end; 


function phidot_ji(t : float) : float; 
var tdj : float: 
begin 
tdj := thetadot_j(t); 
Phidot_j := tdji-(m*wor); 
enc; 


function phi_o(t +: float) : float; 
var tho : float; 
begin 
tho := theta_o(t); 
ohi_o := tho-(m*worst); 
end; 


function phidot_o(t : float) : float; 
var tdo : float; 
begin 
tdo := thetadot_o(t); 
phidot_o := tdo-(m8w_r); 
end; 


function w_i(t =: float) : float; 
var betat, 00+ Pi» PdoOr Pdje Numere denom : float; 
begin 
betat := betal(t); 
O10) == pido e):: 
Pi sssenie eer. 
pdo phidotsoCe jis 
pd4 phidot_i(t); 
numer := betat*betat#*pdjrodo+(pdjrodo)*betat*cos(pj-po); 
denom := l*2*betat*cosl(oj-po)+betat*betat; 
wlio ose m*wore(numer/denom); 
end: 


{SF-} 


mu it 


{SFe) 
function TNTarsetF(t : float) : float: 


var templ : float; 


oegin 

tenpl <= w s(t); 

TNTarsgetF := templ*r*cos(w_ort); 
end: { function TNTersetF } 
{S$F-} 


CSE} 
function TNIargetF_sin(t : float) : float; 


var templ : float: 


begin 

temol s= wilt): 

TNTargetF_sin := templ*r*sin(w_r™t); 
Ne { function TNTarsetF_sin } 
oF} 


procedure Initialize(var LowerLimit : float: 
var UpperLimit =: float; 
var Integral : float; 
var Tolerance : float; 
var Mavintervals : integers; 
var NumIntervals : integer; 
var Error : byte); 


ar a rr enn enn en en eeeren= ) 
{- Output: LowerLimit. UopterLimit. Integral, Tolerance, -) 
{- Maxintervals. Numintervaas» Error -) 
(~ -) 
{- This procedure initializes the acove variables to zero -) 
(a a re treet et eccce= } 


36 


phic sys 


writeln; 
LowerLimit ;= 
UpperLimit : 
Integral ;:2 03 
Tolerance ;:=2= 0 


MaxIntervals := 0; 
Numintervals := 0; 
Error s= 0; 
end: = { procedure Initialize } 


procedure Re_Initializel 
var Intesral : float; 
var NumIntervals ;: integer; 
var Error : byte); 


begin 
Intesral ;:= 0; 
NumIntervals : 
Eamon c=) (0): 
end; { procedure Re_Initialize } 


= 0) 


Procedure get_root_name_of_files(var FileName : strings; 
default_filename : string); 
besin 
writeln: 
FileName := default_filename; 
write('Enter root of file name [°, FileName, °) °); 
readin(FileName): 
if FileName = °" then FileName := default_filename; 
end: { procedure GetOutputFile_Name } 


procedure GetOutputFile_Name(var OutFile : text; 
var FileName : strings; 
default_filename : string); 
var 
Ch : char: 


besin 
repeat 

Chiess) 2 Vi7: 
writeln: 
FileName := default_filenames 
writeC*’Enter file mame [°.» FileNamer °] ° 33 
readin(FileName): 
if FileNanre = ‘'* then FileName := default_filenene; 


assign(OutFile» FileName); 
reset(OutFile): 
if foresult = 0 then { The file already exists. } 
besin 
close(OutFile); 
writelns: 


write('Tnis file already exists. °); 
write('Write over it (Y/N)? °); 

Ch := upcase(ReadKkey); 

writeln(Ch): 

if Ch = ‘¥° then assisn(OutFi le, FileName); 


end: 
ich s vo then 

besin 
rewrite (OQutFile); 
10Check:; 

end: 

Unita CtGhi = 9V¥*) and notl(loerr) i: 
end: { procedure GetOutputFile_ Name } 


procedure GetData(var LowerLimit : float: 
var UpperLimit : float: 
var Tolerance : float: 
var Maxintervals ;: inteser): 


Eee 2 2 2 2 228 2S OSS O22 SC C2 2 SS OS Tee eee ere ez eeeeoeaneeoaen } 
is Output; LowerLimit, UpperLimit, Tolerance. MaxIntervals a 
{- This procedure assisns values to the above variables -} 
{- from keyboard input -} 
(nr rr rn rr nr nr reer en ween cer enn~ } 
procedure GetLimits(var LowerLimit : float; 
var UpperLimit : float): 
a a a ee a me me a ew a er me OOo } 
{- Output: LowerLimit, UpperLimit =) 
i= =}} 
{- This procedure assigns values to the limits of -} 
{- integration from keyboard input a 
besin 
repeat 
repeat 
LowerLimit := 0.03 
1OCheck; 
until not 1O0errs 
weiteln; 
repeat 
UpperLimit := 2™pi/w_or; 
10Check; 
until not 1Oerr;: 
if LowerLimit = UpperLimit then 
besin 
writeln; 
writeln(' The limits of integration must be different.’ ); 
writeln; 
end: 
untii LowerLimit <> UpperLimit: 
end; { procedure GetLinits } 


procedure GetTolerance(var Tolerance : float); 


{- Output: Tolerance -} 


oy 


te -) 
{= This procedure reads in the accepted Tolerance -} 


{= from the keyboard. -} 
Se eee eee noe one ee eee eee ee Beene nee eee } 
const default_tolerance = le~2; 
begin 
writeln; 
repeat 


Tolerance := default_tolerance; 
write(’Tolerance in answer: (> 0): '); 
ReadFloat(Tolerance); 
TOCheck: 
4f Tolerance <= 0 then 
begin 
lOerr := true; 
Tolerance := default_tolerance; 
end; 
Unite net. lOentn: 
end; { procedure GetTolerance } 


procedure GetMaxIntervals(var MaxIntervals : inteser); 


{- Gutput: MaxIntervals - 


{- This procedure reads in the maximum number of - 
{- subintervals to be used in approximating the = 
{- integral. Input 18 from the keyboard. -} 


begin 
writelns 
repeat 
MaxIntervals := 1000; 
write (’Maximum number of subintervals (> 0): °)s3 
ReadInt(MaxIntervals): 
TOCheck: 
4f MaxIntervals <= 0 then 
begin 
ICerr := true; 
MaxIntervals := 1000; 


end: 
untill not JOerr; 
end; { procedure GetMaxiIntervals } 
begin { procedure Ge'Data } 


GetLimits(LowerLimit,. UpperLimit); 
GetTolerance(Tolerance); 
GetMarvIntervals(MaxIntervels); 
{ GetOutputFile(OutFile?); } 

end; { procedure GetDate } 


procedure Results(LowerLimit : float; 
UpperLimit <=: floats 
Tolerance : float; 
MaxIntervals : integer; 
Integral : float; 
NumIntervals :; integer: 
Error : byte}; 


begin 
writeln(OutFile?; 
writeln(OutFile, *Lower LimJt:°:25, LowerLimit:25): 
writeln(QutFile, ‘Upper Limit:°:35,. UpperLimit:25); 
writeln(OutFile. *Tolerance:°:35. Tolerance:25); 
writeln(QutFile, °"Mayimum number of subintervals:’:25, MaxIntervals:5); 
writeln(QOutFile. ‘Numceer of subintervals used:’:35, NumIntervals:5); 
Tee rroOr = cs atihen 
DisplayWarning; 
4¢ Error in [). 2] then 
DisplayError:; 


case Error of 
0 : writeln(OutFile. ‘Integral:':25, Integral}; 


1} : writeln(OutFile. ’Tnhe tolerance must be greater than zero.’ )3 


2: writeln(OutFile, 
"The maximum number of Intervals must be greater than Zero.’ ); 


3: begin 
writeln(QutFile. ‘The integral was not found with ‘'» NumiIntervals, 
’" subintervals.’ ); 
writeln(OutFile, ‘The integral thus far: ‘, Integral); 


end: 
end: { case ) 
end; { procedure Results } 
procedure get_betea_jJ{ver beta_j : float); 
besin 
writeln; 
repeat 


beta_j := beta_j_default; 
write betas jis) 3); 
ReadFloat(beta_J); 
IoCheck:; 
if beta_4j < 0 then 
begin 
IGOerr := true: 
beta_j := beta_jJ_default; 
end: 
until mot JIOerr; 
end; { procedure get_beta_j } 


procedure get_del_theta(var del_theta : float); 
begin 
writelns 
reoeat 
del_theta := del_theta_default; 
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write(*del_thete <: *); 
Resd=loat(del_theta): 
HOSS ck: 
if dei_theta < 0 then 
begin 
lOerr += true: 
del_thete := del_thete_defeult: 
end: S 
until not iOerr;: 
end; { procedure get_del_thete )} 


Procedure get_rho_j (var rho _j : floet): 
begin 
writeln; 
repeat 
rho_j := rho_j_default; 
iG es Cacti Ou 4 mucin is 
ReadFloat(rho_j);: 


MOG eICIY 
aoe. J < 0 then 
begin 
Mei Ro iP iirG 
Fno_j :] rho_}default: 
end; 
Wires al Talelee Niels ele 
hoon start :=s rhold: 
end; { procedure get_del_theta } 
procedure sget_thetec_jJ(ver thetec_j : float); 
besin 
writeln: 
repeet 


BHemeces 2s the tace so default: 
earliest “hietaiGm) (2 (7) 
ReadFlcat(thetac_j):; 

Te check: 

if thetac_jJ < 0 then 


vesgin 
1Oerr := true: 
anetac. J) 22 thetec J default: 
end: 
untii mot JOerr: 
end; { procedure get_del_theta } 
procedure get_del_ wor _revivar del_wir_rev : float): 
begin 
Hip ot Ve Ie 
repeat 


del_w_"_rev := del_w_r_rev_default: 
MiGRee OC am meV sm > )is 
ReacFioat(del_w_r_rev); 
1OCheck: 
fede) wir rey < 0 then 
begin 
NO Siem: Seat Wels 
deil_wior_rev := del_w_or_rev_defeult; 
end: 
mata) not 1Gerr: 
enc: { procedure get_del_thetsa } 


procedure init_pacams; 


begin 
set_beta_j3(€beta_ 4): 
Be. ce, thevatcel. theta): 
Beeman oe CaING. 3) 
get_thetac_jlthetac_J); 
get _del_w_r_revidcel_wior_rev):; 
Wy fe Bee ets Sfalnee 
Cem vame = CC lewecuh evita «pit; 
(SSW GS ecstasy Mies 
feb) BP te Otol RS Tee 
CGE Gaeteeei= tial / WinlGan $15 
pts = te ticle Juan die 
ei is 8 9S telfash. Shige “ayer le os 
Ee Ms ee we Oe lt Ss 
VED eC Mm tmCy ie CO UM Lim ss U lim.t eyelids 
Hi SkeyGa) Qed (ets volt U ee ae 
Voeeee tere < OC de i jsm <= Lilim+et_eyc_i; 
rie BE lal Vina te 3 
syecac | 2= thetac_T init: 


end: 


procedure urite_one_iteration_line(var out_file : tert): 


besin 
writelin(out_file, 

1:4. 
‘ ae minions: Gi Ss 
0 Fs ielels 7s So 
Yo *einiStsEle Wey Baie 
Ui ede 8 SE ca laren sre 
es Aisi ris5 « 
ms fange ©ct+7: 5 « 
TO ASRIGASE) heise 
. as News 7 2S 
Ys, (Sh MI OS NE 

end; 


procedure csj3_seeker(time ;: float; 
var seeker_theta_prev, 
ceever_pni ncrev. 
SCG rer meme tid Scull. 
seerer chji_cur : float); 
Vidi One YC ViG@yCsschn yd © £.0a%s 


ou 


begin 
UrperLimit := time: 
LowerLimit := UpperLimitet_cyc: 


*p := Seeker theta _prev-foca, length; 
yp := seeker_phi_prevefocal_iength: 
ye = seeker_theta_cur*’focal length; 
ye 3:= seeker_phi_cur*focal_lensth; 

xd 32 ype-~C? 

yd :2 yp-yc; 

rho 3:2 sart({vd«xd)+lyd-yd)); 
thetac_T := atanedlyd. xd}; 


edeptive Gauss _Quadrature(LowerLimit. UpperLimit. Tolerance. MaxIntervals, 
Intecrale NumIntervals» Errors. oTNTargetF }; 

“x td .- J[ntesralstoiem): 

Re_Inijtialize(Integral. Nurintervals. Error): 

Adaptive_Gauss_Quacrature(LowerLimit. UpperLimit.s Tolerance. Maxintervels, 
Integrals NumIntervals. Error. aTNTarsetF_sin); 

y_tj :2# Integral/s(piem); 


{ renge_tot:=sart(r_ti*«<_ ti+y_ti*y_ti); 
theta_tgt:satanily_tj>»~_t3); 
oel_range:=range tat: 
del_x:ecel_ranoecos(theta_tgt): 
del_y:2del_ranse“sin(theta_t9t);)} 


we te 


xC 22 yvp-del_x: 

ye r= ye-cel_y: 

seeker_theta_cur 3 

seerer_phi_cur :s 

seeker_t"eta_rrev 

seeker_phi_prev : 
enc: 


yes focel_lenath: 
cefocai_lensgth: 

= seeker_thete_cur: 
seeker_phi_cur: 


procedure write _table_header(var out_file : text): 
begis 

wreiteln(out_file?}: 

writein(out file. 


’ 
,’ 
’ 
' 
’ y_¢3? 
’ 
’ 
’ 
’ 


end: 


procedure write_end_iterations_line(ver out_file =: text); 


begin 
writelnCout_file. 

1:6. 
Cee ili) aio Sc 
Po Metaleig Sic 
eae tid Clalncen cron 
Yet elo on oe 

end: 


Frocedure write _file_Neader: 

begin 
writeln( 
writelnt 
7 ey Un A 
wristelnt 


AON 
109 0 
te 
Tt 
* 
733 
sof =e 
soe, IN 
oan and wie: 
ta EN en 
Le et | ee) | 
we we 


Oe 


=e 


oe ae a ee 


wrdtein(esj 
WElewedin Ces sie 
writeln(css_ 
writelnat esd. 
Crivelmvas 
writeln(csi_file. ted 
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